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Abstract 

We  numerically  examine  the  effectiveness  of  thermoacoustic  actuators  that  locally 
introduce  high-intensity  acoustic  waves  for  active  flow  control.  The  use  of  near-field 
acoustic  excitation  in  flow  control  applications  has  not  been  explored  deeply  in  the  past 
due  to  the  absence  of  surface-compliant  actuators  that  are  capable  of  delivering  high- 
intensity  acoustic  waves.  In  this  investigation,  we  perform  LES  of  flow  control  with 
acoustic  waves  motivated  by  the  recent  development  in  manufacturing  graphene/carbon 
nanotube-based  surface  compliant  loud  speakers.  The  interaction  of  the  acoustic  waves 
and  turbulent  flow  structures  is  analyzed  with  high-fidelity  LES  for  compressible  flow 
over  a  wall-mounted  hump  at  a  Reynolds  number  of  0.5  x  106  and  Mach  number 
of  0.25.  In  the  computations,  we  consider  the  use  of  high-frequency  actuation  at 
Helmholtz  number  of  3  based  on  performance  characteristics  of  the  graphene-based 
thermoacoustic  actuators.  We  observe  that  the  actuation  is  able  to  introduce  small- 
scale  perturbation  to  the  shear  layer  in  the  separated  flow  and  delay  the  formation  of 
large-scale  spanwise  vortices.  This  hence  elongates  the  recirculation  zone  and  shifts 
the  low-pressure  region  downstream  of  the  hump.  As  a  consequence,  the  drag  on  the 
hump  is  reduced  by  approximately  5%  for  the  three-dimensional  calculation  and  up  to 
10%  in  the  two-dimensional  flow.  The  mechanism  of  drag  reduction  is  different  from 
the  one  using  synthetic  jets,  which  relies  on  reducing  the  size  of  the  recirculation  zone. 
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1  Major  Accomplishments 


1.  Conducted  large-eddy  simulations  (LES)  of  three-dimensional  turbulent  flow  over  a 
wall-mounted  hump  (NASA  validation  workshop  problem)  at  Re^  =  0.5  x  106  and 
Mqo  =  0.25  and  validated  the  results  with  other  experimental  and  computational 
studies.  This  computational  setup  serves  as  a  numerical  test  bed  for  various  flow 
control  parameter  studies. 

2.  Formulated  a  thermoacoustic  actuator  wall  model  that  be  used  in  LES  or  DNS  to 
predict  actuator  performance  in  flow  control  applications. 

3.  Performed  two  and  three-dimensional  flow  control  simulations  with  thermoacoustic 
actuation  modeled  through  a  wall  boundary  condition.  Modification  of  the  turbulent 
flow  field  and  reduction  in  drag  on  the  hump  are  observed. 

4.  Analyzed  the  interaction  between  the  turbulent  separated  flow  and  acoustic  waves 
responsible  for  modifying  the  wake  downstream  of  the  hump. 

2  Introduction 

Active  flow  control  has  been  applied  to  various  aerodynamic  systems  for  performance  en¬ 
hancement,  such  as  lift  increase,  drag  reduction  and,  mixing  enhancement  [1].  The  main 
challenge  of  active  flow  control  is  to  determine  how  and  where  one  can  introduce  the  control 
input  to  modify  the  flow  field  favorably  with  minimal  energy  input.  In  past  applications, 
active  flow  control  has  been  performed  by  injecting  mass,  momentum,  or  vorticity  (or  the 
combination  thereof)  in  either  a  steady  or  unsteady  fashion  using  various  types  of  devices 
[2],  such  as  the  synthetic  jets  [3],  dielectric  barrier  discharge  (DBD)  plasma  actuators  [4], 
or  combustion  powered  actuators  [5].  Compared  to  passive  flow  control  techniques,  such 
as  vortex  generators,  the  use  of  active  flow  control  has  not  been  widespread  for  general 
aerodynamic  enhancements  (despite  the  adaptive  capability  of  active  flow  control).  There 
are  two  main  factors  for  restricting  the  implementation  of  active  flow  control.  First  is  the 
bulkiness  due  to  the  weight  and  size  of  the  actuator  device  which  limits  the  location  it  can 
be  embedded.  Second  is  the  amount  of  power  required  to  deliver  the  control  perturbation 
into  the  flow  field,  which  is  fairly  large  because  of  the  limited  device  efficiency. 

To  address  the  above  two  limitations,  we  consider  the  use  of  a  novel  speaker  made 
of  carbon  nanotubes/graphene  as  a  flow  control  actuator.  A  membrane  of  carbon  nan¬ 
otubes/graphene  laid  on  a  backing  material  such  as  paper  under  AC  power  has  been  shown 
to  produce  large-amplitude  acoustic  waves  [6] .  The  input  AC  current  sinusoidally  heats  this 
device  due  to  joule  heating  and  creates  surface  pressure  disturbances  that  are  emitted  as 
acoustic  waves  to  the  surrounding  fluid  [7,  8].  One  of  the  attractive  characteristics  of  this 
carbon-based  sound-emitting  device  is  its  sheet-like  arrangement  that  can  be  applied  to  es¬ 
sentially  any  surface,  without  occupying  any  internal  space  nor  adding  significant  weight, 
which  is  similar  to  how  DBD  plasma  actuators  can  be  installed. 
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The  sound  generation  mechanism,  known  as  thermophone,  was  discovered  with  platinum 
much  earlier  (1917)  [9]  but  has  not  been  widely  utilized  due  to  the  relatively  large  power  re¬ 
quired.  Compared  to  conventional  metals,  the  heat  capacity  value  for  carbon-based  material 
(carbon  nanotubes/graphene)  in  consideration  here  is  at  least  2  orders  of  magnitude  smaller. 
Since  the  output  acoustic  power  delivered  to  the  surrounding  flow  field  is  related  inversely 
to  the  material  heat  capacity  C  (i.e.,  P0,,tp„t  oc  C'_1Pinput),  the  use  of  carbon-based  material 
dramatically  enhances  the  device  efficiency.  According  to  the  measurements  reported  in  Tian 
et  al.  [6]  and  Lin  et  al.  [10] ,  the  carbon-based  devices  are  able  to  deliver  constant  noise  levels 
of  ~  85dB  over  a  broad  range  of  frequency,  10-50  kHz,  which  should  be  of  interest  for  many 
flow  control  applications  [11], 

Past  studies  have  shown  that  inherent  acoustic  feedback  mechanisms  in  flow  over  a  cavity 
[12]  or  impinging  jet  [13]  yield  significantly  different  flow  behavior  compared  to  cases  without 
the  presence  of  such  feedback  path.  There  are  also  reports  of  acoustic  perturbation  being  able 
to  modify  the  flow  over  airfoils  [14]  or  stabilize  combustion  instabilities  [15]  from  external 
acoustic  input  (not  attached  to  the  body  surface).  We  however  re-emphasize  that  past 
acoustic  actuators  could  not  be  readily  mounted  on  the  surface  of  bodies  due  to  their  size 
(including  synthetic  jets  that  relies  mostly  on  momentum  injection).  The  present  study  will 
examine  how  strong  acoustic  perturbation  introduced  locally  can  alter  the  flow  field  based 
on  a  surface-mounted  thermoacoustic  actuator. 

The  objective  of  the  present  study  is  to  gain  insights  into  how  thermoacoustic  actua¬ 
tors  based  on  carbon  nanotubes/graphene  can  alter  turbulent  separated  flows.  A  boundary 
condition  model  [16]  that  captures  the  generation  of  thermoacoustic  waves  is  implemented 
in  the  LES  for  flow  over  a  wall-mounted  hump  to  examine  the  influence  of  thermoacoustic 
perturbation  on  reducing  drag  and  shortening  the  recirculation  bubble.  The  hump  geometry 
is  chosen  from  the  NASA  Validation  Workshop  [17,  18,  19,  20],  shown  in  Figure  1.  The  use 
of  this  hump  geometry  facilitates  the  assessment  of  the  thermoacoustic  actuator  performance 
against  other  types  of  actuators.  This  specific  problem  is  chosen  to  assess  the  performance 
of  high-intensity  acoustic  waves  for  turbulent  flow  modification.  The  Reynolds  number  and 
Mach  numbers  are  chosen  to  be  Re  =  0.5  x  106  and  =  0.25,  respectively. 


Figure  1:  Geometry  of  the  NASA  workshop  hump  model. 


Below  we  describe  computational  approach  taken  in  this  study  and  present  the  validation 
cases  in  Section  3.  We  then  describe  the  thermoacoustic  actuator  in  Section  4.  Flow  control 
findings  from  the  performed  LES  for  two  and  three-dimensional  flow  are  presented.  It  is 
found  that  the  thermoacoustic  actuator  can  alter  the  turbulent  flow  field  and  achieve  drag 
reduction  on  the  wall-mounted  hump  in  Section  5.  At  last  in  Section  6,  we  provide  a 
summary  of  the  present  study  and  comments  on  possible  future  tasks  to  further  understand 
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the  mechanism  of  active  flow  control  using  locally  added  high-intensity  acoustic  waves. 

3  Computational  Approach 

Numerical  simulations  of  flow  over  the  wall-mounted  hump  are  performed  with  CharLES 
[21,  22,  23]  to  accurately  predict  the  turbulent  flow  field  as  well  as  the  acoustic  field.  The 
accuracy  of  the  compressible  flow  solver  CharLES  is  formally  second  order  in  space  and 
third  order  in  time.  The  flow  over  the  hump  is  examined  using  LES  with  the  Vreman  model, 
since  the  small-scale  vortical  structures  generated  at  the  considered  Reynolds  number  can  be 
smaller  than  the  grid  resolution  near  the  wall.  The  Reynolds  number  and  Mach  number  of 
Re  =  PooUoqc/ Hoo  =  0.5  x  106  and  M0 0  =  U^/a <*,  =  0.25  are  chosen  to  match  representative 
cases  considered  in  the  NASA  Validation  Workshop  [18,  19,  20]. 

The  computational  domain  size  is  chosen  to  match  the  experimental  setup  given  by  the 
workshop  guidelines  [17]  and  the  study  conducted  by  Greenblatt  et  al.  [24],  The  streamwise 
extent  of  the  computational  domain  is  x/c  €  [—2.4, 4]  with  the  hump  positioned  from  x/c  =  0 
to  1,  as  illustrated  in  Figure  2.  The  vertical  extent  of  the  domain  is  y/c  =  0.909  with  the 
top  boundary  lowered  to  contract  the  flow  above  the  hump.  This  contraction  is  introduced 
to  capture  the  effects  of  the  side  splitter  plates  [17].  The  spanwise  direction  is  taken  to  be 
periodic  with  an  extent  of  z/c  €  [0,  0.2],  The  spatial  domain  is  discretized  with  a  structured 
grid  with  a  near-wall  grid  resolution  of  Ax+  =  utAx/is  =  25,  Ay+  =  0.25,  Az+  =  15  around 
the  humD.  The  overall  crid  size  for  the  baseline  ( uncontrolled )  simulation  is  48  x  106. 


Figure  2:  Geometry  of  the  computational  domain  (xy-plane).  The  spanwise  (z)  extent  is  0.2  for  three- 
dimensional  computations.  Turbulent  inflow  condition  is  specified  at  the  inlet.  Up-winding  and  the  sponge 
zone  are  utilized  to  allow  waves  to  leave  the  domain  without  numerical  reflections. 

This  particular  turbulent  flow  simulation  is  known  to  be  sensitive  to  the  upstream  condi¬ 
tion.  The  inlet  boundary  condition  is  set  to  match  the  experimental  turbulent  velocity  pro¬ 
file  reported  in  the  workshop  guidelines  [17]  and  the  experiments  performed  by  Greenblatt 
et  al.  [24] .  The  inlet  boundary  layer  profile  is  specified  to  match  experimental  measurement 
profile.  We  also  ensure  that  the  value  of  turbulent  intensity  is  matched  at  the  inlet  to  allow 
for  turbulence  to  develop  before  reaching  the  leading  edge  of  the  hump.  Turbulent  fluctu¬ 
ation  is  added  to  the  inlet  profile  with  random  Fourier  modes  [25,  26].  No-slip  boundary 
condition  is  applied  along  the  wall  boundary.  Along  the  outlet  and  the  top  of  the  bound¬ 
ary,  we  utilize  the  sponge  zone  to  damp  out  any  outgoing  acoustic  waves  or  exiting  vortical 
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structures  [27].  A  first  order  up- winding  scheme  is  applied  away  from  the  hump  to  damp 
out  the  exiting  waves  to  ensure  that  waves  do  not  artificially  reflect  back  into  the  near  field. 
The  top  boundary  also  utilizes  the  no-stress  boundary  condition. 

3.1  Validation 

To  ensure  that  the  present  calculation  accurately  captures  the  flow  physics,  we  compare 
the  baseline  separated  flow  over  the  hump  with  the  experimental  measurements  of  Seifert 
and  Pack  [18]  ( =  0.25,  Re^  =  2.4  to  26  x  106),  Greenblatt  et  al.  [24]  (. M, ^  =  0.1, 
Re oo  =  1  x  106)  and  the  numerical  results  of  Franck  and  Colonius  [20]  =  0.25,  Re ^  = 

0.6  x  106).  The  present  results  are  compared  with  the  aforementioned  results  in  Figures 
3,  4,  and  5.  The  pressure  coefficient  cp  =  {jp  —  Poo) / {\PooU^)  is  compared  in  Figure  3.  It 
can  be  observed  that  the  present  baseline  result  agrees  very  well  with  other  experimental 
measurements  throughout  the  domain.  Our  result  also  agrees  reasonably  well  with  numerical 
calculation  by  Franck  and  Colonius  [20]  but  we  note  that  the  present  study  uses  a  larger  grid 
size  than  theirs,  which  might  have  contributed  to  the  slight  difference. 


x/c 

Figure  3:  Comparison  of  pressure  distribution  over  the  hump. 

The  corresponding  time-average  velocity  profiles  are  presented  in  Figures  4  and  5  and 
agree  well  with  the  PIV  measurements  of  Greenblatt  et  al.  [24]  and  computational  findings 
of  Franck  and  Colonius  [20].  The  Reynolds  stress  distribution  is  also  found  to  match  the 
experimental  measurements.  From  these  results,  it  has  been  observed  that  there  are  many 
common  traits  in  the  separated  flow  at  these  Reynolds  numbers  {Re ^  =  0.5  x  106  and  I06) 
and  Mach  numbers  (M0 0  =  0.1  and  0.25).  There  is  some  elongation  of  the  recirculation 
region  downstream  of  the  hump  for  higher  Mach  number  ~  0.6  as  reported  in  Franck  and 
Colonius  [20]. 

For  the  chosen  grid  resolution  and  domain  setup,  the  baseline  results  are  found  to  be  in 
agreement  with  those  reported  in  past  studies.  The  validated  baseline  simulation  provides 
us  with  the  test  bed  to  explore  the  use  of  thermoacoustic  flow  control  on  the  wall-mounted 
hump  problem. 
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Figure  4:  Comparison  of  time-average  velocity  and  Reynolds  stress  for  the  baseline  three-dimensional  flow. 
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Figure  6:  Uncertainty  quantification  of  the  model  P  =  S  (a\T  +  0L2 Q  +  a3  7jr)  ar|d  its  (a)  comparison  with 
the  experimental  measurements  of  Tian  et  al.  [6].  The  shaded  region  constitutes  the  propagated  uncertainty 
from  the  probability  density  profile.  The  correlations  (b)  between  the  coefficients  and  their  probability 
distribution  (c)  are  also  revealed  with  MCMC. 


4  Thermoacoustic  Actuator  Model 

The  generation  of  acoustic  waves  from  graphene/carbon-nanotube  based  actuators  is  mod¬ 
eled  by  a  thermal  boundary  condition  in  the  present  LES.  In  the  early  work  of  Arnold  and 
Crandall  [9],  the  input  electrical  power  P  is  related  to  the  temperature  T  of  the  actuator 
and  its  temporal  derivative  P  =  S  (b\T  +  627^7  ),  where  the  coefficients  are  dependent  on 
the  material  and  configuration  and  S  is  the  surface  area  of  the  actuator.  The  first  term 
represents  convective  heat  transfer  and  the  second  term  captures  the  effect  of  capacitance. 
In  addition  to  the  above  formulation,  there  are  models  [28]  that  include  an  extra  heat  flux 
term,  P  =  S  {a{T  +  a2q  +  a3^) . 

We  have  an  ongoing  theoretical  and  numerical  investigation  to  quantify  the  significant 
of  each  term  in  the  above  model.  The  objective  of  this  ongoing  effort  is  to  quantify  the 
acoustic  contribution  from  each  term  to  accurately  predict  the  performance  of  thermoacoustic 
actuators.  This  question  is  approached  with  the  MCMC  (Markov  Chain  Monte  Carlo) 
method  [29]  to  determine  whether  there  is  a  correlation  between  the  coefficients  and  how 
sensitive  the  choice  of  coefficients  are  since  they  do  not  appear  to  be  directly  measurable. 
Preliminary  computations  are  shown  in  Figure  6  that  compare  the  experimental  measurement 
of  Tian  et  al.  [6]  and  the  present  MCMC  results  with  uncertainty  quantification. 

For  the  purpose  of  the  current  flow  control  simulations,  we  have  found  that  the  use 
of  sinusoidal  heat  flux  g"all  =  ecos(aA)  (where  e  <  1)  to  represent  the  actuator  predicts 
the  characteristics  of  the  emitted  acoustic  waves  reasonably  well,  as  shown  in  Figure  7.  A 
companion  study  [16]  examining  this  actuator  model  has  extended  the  acoustic  piston  model 
[30].  We  have  found  that  the  thermal  wavelength  is  very  small  compared  to  the  acoustic 
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Figure  7:  Simulated  spherical  thermoacoustic  waves  generated  from  a  sinusoidal  heat  flux  input  on  a  flat 
plate.  DNS  results  (left)  and  comparison  of  directivity  with  the  FSU’s  analytical  model  (right). 


wavelength  and  the  thermal  wave  decays  exponentially  with  the  distance  from  the  wall.  We 
emphasize  that  the  current  flow  control  experiments  are  mainly  concerned  with  how  the 
acoustic  waves  interact  with  the  turbulent  separated  flow  and  does  not  focus  deeply  into 
how  the  waves  itself  is  created.  Hence,  we  use  the  above  heat  flux  boundary  condition  in  the 
active  flow  control  simulations  shown  below.  We  note  that  the  ongoing  study  using  MCMC 
will  play  an  important  role  in  determining  the  mechanism  of  how  thermoacoustic  actuators 
generate  high-intensity  waves  and  will  guide  the  optimal  design  of  these  actuators  down  the 
road. 

5  Active  Flow  Control  with  Acoustic  Waves 

In  flow  control  simulations,  we  position  the  thermoacoustic  actuator  (boundary  condition) 
near  the  separation  point  as  a  spanwise  strip  with  a  width  a  of  cr/c  =  4.4  x  10-3.  The  width 
is  selected  to  be  the  same  as  the  width  of  the  synthetic  jet  actuator  slot  used  in  the  NASA 
validation  study.  The  position  of  the  actuator  is  also  chosen  to  be  at  the  same  location, 
which  is  x/c  =  0.6563.  In  future  analysis,  the  importance  of  the  actuator  size  and  location 
will  be  examined.  For  the  controlled  cases,  we  must  resolve  the  interaction  between  the 
high-frequency  acoustic  waves  emitted  by  the  actuator  and  the  surrounding  turbulent  flow. 
This  required  us  to  further  resolve  the  original  grid  and  increased  the  overall  grid  size  to  be 
94.4  x  106  for  the  three-dimensional  controlled  case. 

The  primary  frequency  used  in  this  study  was  based  on  a  Helmholtz  number  ( ka  =  3)  [6] . 
For  the  experimental  study  the  range  of  desirable  frequencies  for  outputting  high-intensity 
acoustic  waves  are  above  10  kHz  with  the  graphene-based  actuator  of  size  ~  lcm.  In  this 
study,  we  normalize  the  frequency  based  on  the  width  of  the  actuator  ((9(10-3))  such  as 
fa  =  f  c/ cIqq  =  107,  where  /  is  the  dimensional  frequency.  For  the  two-dimensional  case,  we 
also  consider  the  higher  harmonics  of  2  f0  and  IO/q. 

The  validity  of  the  heat  flux  boundary  condition  is  examined  on  the  hump  in  a  two- 
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Figure  8:  2D  simulated  thermoacoustic  waves  generated  from  a  sinusoidal  heat  flux  input  on  the  hump. 
LES  pressure  contour  (left)  and  its  spectrum  at  5w  away  from  the  actuator.  SPL  is  112  dB  at  this  location. 


dimensional  domain  with  a  quiescent  condition  using  an  actuation  frequency  of  u>  =  2irfa  = 
27t(107)  and  an  amplitude  of  eq  =  0.00153,  as  shown  in  Figure  8  (left).  The  pressure  is 
also  monitored  at  the  location  5cr  away  normally  from  the  hump  surface  and  its  spectrum 
is  shown  in  Figure  8  (right).  The  measured  SPL  at  this  location  is  ~  110  dB,  which  is 
larger  than  ps  85  dB  as  reported  by  Tian  et  al.  [6],  who  took  measurements  from  a  square 
area  of  actuator  rather  than  a  stripe  arrangement.  We  use  this  large  amplitude  of  HOdB  to 
ensure  we  observe  difference  in  the  flow  with  control.  Follow-up  investigation  will  consider 
lower  actuation  amplitude  and  examine  its  influence  on  the  control  effectiveness.  We  note 
in  passing  that  the  drag  and  lift  (coefficients)  directly  generated  by  actuation  in  quiescent 
flow  are  of  order  10— 6 ,  which  are  practically  negligible. 

We  present  findings  from  the  two  and  three-dimensional  LES  with  thermoacoustic  bound¬ 
ary  condition  applied  for  active  flow  control.  Below  we  contrast  the  difference  between  the 
uncontrolled  and  controlled  turbulent  flow  fields  and  discuss  how  the  exerted  aerodynamic 
force  is  altered.  The  three-dimensional  LES  requires  significant  amount  of  computational 
resource  for  the  Reynolds  number  of  0.5  x  106  since  the  grid  sizes  are  large  (48  and  94  mil¬ 
lion  points  for  the  baseline  and  controlled  cases,  respectively).  While  we  are  able  to  perform 
only  a  limited  number  of  three-dimensional  computations  due  to  the  large  computational  re¬ 
source  required,  we  have  also  conducted  two-dimensional  companion  simulations  to  examine 
if  acoustic  flow  control  can  alter  the  flow  field  with  only  two-dimensional  instabilities.  The 
acoustically  controlled  flow  shows  favorable  results  in  reducing  the  drag  on  the  wall-mounted 
hump  for  both  two  and  three-dimensional  cases. 

5.1  Two-Dimensional  Flow 

First,  we  consider  flow  control  with  high- intensity  acoustic  waves  for  two-dimensional  flow. 
The  computational  domain  is  chosen  to  be  the  same  for  three-dimensional  simulations  but 
with  only  a  single  cell  in  the  spanwise  direction  (no  spanwise  variation).  Simulations  with 
three  actuation  frequencies  of  fa  =  107,  214, 1070  are  performed,  in  addition  to  the  baseline. 
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Figure  9:  Instantaneous  vorticity  contour  of  the  two-dimensional  controlled  flow  for  different  actuation 
frequencies. 


Figure  10:  Numerical  Schlieren  of  the  two-dimensional  controlled  flow  for  different  actuation  frequencies. 


The  actuation  frequency  in  terms  of  the  Helmholtz  number  is  approximately  3  for  fa  =  107, 
which  is  based  on  the  cutoff  frequency  achieved  in  the  experimental  settings  of  Tian  et  al.  [6]. 

Due  to  the  absence  of  spanwise  variation,  two-dimensional  simulation  reveals  larger  span- 
wise  vortices  forming  and  shedding  from  the  wall-mounted  hump.  The  spanwise  vorticity 
contours  are  shown  in  Fig.  9  for  both  the  uncontrolled  and  controlled  cases.  With  high- 
frequency  acoustic  perturbations  added,  the  spanwise  vortices  appear  to  become  smaller 
compared  to  the  baseline  case,  which  is  also  evident  from  the  Schlieren  images  in  Fig.  10.  In 
the  Schlieren  images,  the  emitted  acoustic  waves  can  be  seen  for  fa  =  107  and  214  while  the 
waves  become  hardly  visible  in  the  figure  for  fa  =  1070  for  the  highest  actuation  frequency. 
The  acoustic  perturbations  appear  to  create  small-scale  structures  that  trigger  the  roll  up  of 
smaller  vortices  behind  the  hump.  We  also  observe  the  generation  of  smaller  opposite  sign 
vortices  near  the  wall. 

The  corresponding  time-averaged  vorticity  plots  are  shown  in  Fig.  11,  where  we  observe 
that  spanwise  vorticity  is  transported  (convected)  further  downstream  with  actuation.  We 
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x/c  =  0.8 


Figure  12:  Slice  of  time-averaged  spanwise  vorticity  (—  u2)  at  x/c  =  0.8  for  two-dimensional  flow. 


show  in  Fig.  12  the  time-average  vorticity  profile  at  x/c  =  0.8  to  illustrate  that  the  vortex 
sheet  emanating  from  the  separation  point  becomes  thinner  with  control.  This  allows  for  the 
sheet  to  convect  downstream  with  delayed  growth  of  large-scale  instability  and  the  circulation 
zone  to  elongate  in  the  streamwise  direction,  as  presented  in  Fig.  13.  Also  shown  in  the  same 
figure  are  the  pressure  contours.  The  low-pressure  zone  behind  the  hump  moves  backward 
and  the  recirculation  region  elongates,  in  particular  for  the  controlled  cases  of  fa  =  107  and 
fa  =  1070. 

The  elongation  of  the  recirculation  region  correlates  to  the  cases  with  the  largest  reduction 
in  drag,  as  seen  in  Table  1.  Since  the  low-pressure  core  is  pushed  further  downstream,  the 
drag  force  on  the  hump  is  decreased.  The  pressure  coefficient  over  the  hump  is  plotted  in 
Figure  14.  We  note  that  upstream  of  the  separation  point  there  appears  to  be  only  minor 
differences  between  the  baseline  and  the  actuation  cases.  This  tells  us  that  the  drag  reduction 
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Figure  13:  Time-average  streamlines  and  pressure  for  the  two-dimensional  controlled  flow  for  different 
actuation  frequencies. 


Cases 

<? 

o 

1 

to 

A  CD 

Baseline 

5.56 

— 

fa 

=  107 

5.01 

-9.9% 

fa 
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-0.9% 

fa  : 
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4.93 

-11.4% 

Table  1:  Drag  reduction  with  acoustic  flow  control  for  two-dimensional  flow. 


comes  from  the  modifications  made  in  the  wake.  The  acoustic  flow  control  is  able  to  reduce 
the  drag  on  the  hump  by  up  to  11.4%  for  two-dimensional  flow.  The  wake  elongation  used  for 
drag  reduction  is  of  an  opposite  trend  from  how  synthetic  jet  controls  the  flow  over  the  hump 
in  past  studies.  With  synthetic  jets,  recirculation  zones  are  reduced  in  size  for  recovering 
pressure  losses  to  achieve  drag  reduction. 

Spectrum  of  lift  is  also  shown  in  Figure  15,  where  the  decease  of  dominant  peak  ampli¬ 
tudes  for  the  controlled  case  is  observed.  The  actuation  frequencies  are  also  marked  (see 
logarithmic  plot  on  right)  and  can  be  seen  to  be  of  very  limited  spectral  power.  The  use 
of  acoustic  perturbation  does  not  ‘directly’  shift  the  flow,  but  appears  to  trigger  instability 
differently  to  the  controlled  flow  that  alters  the  global  flow  field.  This  approach  shares  sim¬ 
ilarity  with  the  flow  control  on  a  jet  using  high-frequency  piezoelectric  actuator  performed 
by  Wiltse  and  Glezer  [11]. 

5.2  Three-Dimensional  Flow 

Next,  we  consider  the  effect  of  high- intensity  acoustic  waves  on  three-dimensional  turbulent 
flow  over  the  hump.  Due  to  the  large  computational  time  required  to  perform  the  LES 
calculation  on  the  DoD  cluster,  we  have  limited  our  controlled  simulation  to  be  for  a  single 
case  with  fa  =  107  and  amplitude  of  llOdB  (SPL). 
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Figure  14:  Time-averaged  pressure  coefficient  over  the  hump  for  two-dimensional  controlled  flow. 
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Figure  15:  Lift  spectra  with  flow  control  of  two-dimensional  flow  presented  on  linear  and  logarithmic  scales. 
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Figure  16:  Comparison  of  the  baseline  and  acoustically  controlled  (fa  =  107)  instantaneous  flow  fields 
downstream  of  the  wall-mounted  hump.  Three-dimensional  structures  are  visualized  by  the  Q-criterion 
colored  with  streamwise  velocity.  Note  that  the  turbulent  structures  are  smaller  with  elongated  wake  for  the 
controlled  flow  using  high-intensity  acoustic  waves. 


The  instantaneous  three-dimensional  flow  fields  from  the  baseline  and  controlled  cases 
are  shown  in  Fig.  16.  The  controlled  simulation  required  a  finer  grid  resolution  around  the 
actuator  location  to  resolve  the  interaction  between  the  acoustic  waves  and  the  turbulent 
structures.  The  overall  grid  size  for  the  controlled  flow  is  94  million  points.  We  believe 
both  the  uncontrolled  and  controlled  cases  are  resolved  well  since  the  baseline  case  has  been 
well  validated  as  discussed  earlier  and  the  vortical  wake  structures  for  the  controlled  flow 
are  well-captured  as  it  can  be  seen  in  the  numerical  Schlieren  in  Fig.  17.  We  notice  that 
the  turbulence  flow  structures  become  finer  and  the  wake  is  elongated  with  acoustic  control. 
These  observations  are  aligned  with  those  made  for  two-dimensional  flow  discussed  above. 
The  three-dimensional  vortical  structures  are  of  course  much  finer  due  to  the  variations 
allowed  in  the  spanwise  direction. 

In  Fig.  16,  the  large  blue  region  (low  velocity  region)  is  prominent  in  the  controlled  case 
and  is  obvious  that  separated  region  extends  further  downstream  of  the  hump.  This  growth 
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Figure  17:  Numerical  Schlieren  of  the  three-dimensional  controlled  flow  with  an  actuation  frequency  of 
fa  =  107. 

for  three-dimensional  flow  can  be  seen  clearly  from  Fig.  18.  The  recirculation  zone  for  the 
controlled  case  is  found  to  elongate  by  approximately  15%  compared  to  the  baseline  case. 
This  elongation  also  shifts  downstream  the  low-pressure  region  associated  with  the  wake  and 
leads  to  drag  reduction  as  compared  in  Table  2. 
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Table  2:  Drag  reduction  with  acoustic  flow  control  for  three-dimensional  flow. 

Let  us  further  analyze  the  shear  layer  with  the  time-average  profile  of  spanwise  vorticity 
using  Figs.  19  and  20.  It  can  be  noted  that  the  vortex  sheet  is  thinner  and  convects  further 
downstream  with  the  use  of  actuation.  The  comparison  of  the  vorticity  distribution  shows 
less  spreading  of  vorticity  for  the  controlled  flow  in  Fig.  20.  Thinning  of  the  vorticity  profile 
at  this  location  supports  the  argument  that  the  sheet  needs  to  advect  further  downstream 
to  spread  and  reattach  to  the  wall,  which  in  turn  elongates  the  recirculation  zone.  Since 
the  vorticity  profile  upstream  of  the  separation  point  cannot  change  much,  the  control  effort 
with  acoustics  is  aimed  at  modifying  the  spreading  rate  of  spanwise  vorticity  to  change  the 
wake  profile. 

The  addition  of  the  high-frequency  pressure  oscillations  to  the  flow  creates  small-scale 
structures.  These  small-scale  structures  appear  to  suppress  the  production  of  large  vortices 
that  exist  in  the  uncontrolled  flow.  The  instantaneous  pressure  and  vorticity  contours  in 
Fig.  21  illustrate  this  observation.  In  fact,  the  turbulent  energy  in  the  shear  layer  appears 
to  decrease  with  control  as  shown  by  the  reduction  of  Reynolds  stress  in  Fig.  22.  From  this 
figure,  we  see  that  the  magnitude  of  fluctuations  in  the  shear  layer  decrease  with  the  added 
acoustic  waves.  For  the  unactuated  case,  the  inherent  shear  layer  instability  and  the  vortex 
sheet  roll  up  cause  large-scale  fluctuation  in  the  wake.  Conversely,  for  the  actuated  case 
the  fluctuations  near  the  separation  point  are  decreased  and  the  shear  layer  breaks  down 
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Figure  18:  Time-averaged  streamlines  for  baseline  (top)  and  the  actuated  case  (bottom)  for  three- 
dimensional  flow. 
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Figure  19:  Time-averaged  spanwise  vorticity  (cjz)  for  baseline  (top)  and  actuated  (bottom)  three- 
dimensional  flow. 
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Figure  20:  Slice  of  time-averaged  spanwise  vorticity  (cjz)  for  three-dimensional  flow  at  x/c  =  0.8. 


Baseline  Controlled 


Figure  21:  Instantaneous  three-dimensional  flow  field  images  for  the  baseline  and  actuated  cases.  Pressure 
contours  on  the  top  and  spanwise  vorticity  contours  on  the  bottom. 


further  downstream.  It  appears  that  the  interaction  of  the  high-intensity  acoustic  waves  and 
the  shear  layer  instability  plays  a  large  role  in  how  small-scale  spanwise  vortices  develop 
and  affect  the  formation  of  larger-scale  vortical  structures.  The  reduction  in  the  size  of  the 
vortex  sheet  roll-up  appears  to  result  in  elongating  the  recirculation  zone.  Furthermore,  the 
observations  made  for  three-dimensional  flows  are  in  agreement  with  two-dimensional  flows. 
This  suggests  that  the  main  mechanism  for  acoustic  control  with  the  current  actuator  setup 
is  two-dimensional  in  nature. 


6  Summary 

In  this  investigation,  we  have  numerically  examined  the  use  of  high-intensity  acoustic  waves 
for  the  control  of  separated  flow  over  a  wall-mounted  hump.  The  computations  were  per¬ 
formed  for  a  chord-based  Reynolds  number  of  0.5  x  106  and  Mach  number  of  0.25  using 
large-eddy  simulations.  The  uncontrolled  baseline  calculation  was  validated  with  experi- 
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Figure  22:  Reynolds  stress  at  the  rear  of  the  hump  for  three-dimensional  flow.  Shown  are  the  baseline  (top) 
and  actuated  cases  (bottom). 


mental  and  computational  results  available  in  the  literature  and  was  observed  to  be  in  good 
agreement.  In  order  to  accurately  predict  the  baseline  turbulent  flow  field,  a  grid  size  of  48 
million  was  required,  which  called  for  the  use  of  a  large- size  computer  cluster  made  available 
through  the  HPCMP  program  at  the  Department  of  Defense. 

The  generation  of  thermoacoustic  waves  was  analyzed  to  determine  an  appropriate  bound¬ 
ary  condition  for  LES.  A  sinusoidal  heat  flux  condition  was  found  to  be  sufficient  in  the 
present  study  to  be  utilized  as  a  wall  boundary  condition  to  create  acoustic  waves  having 
sound  pressure  level  of  HOdB,  measured  five  slot  widths  away  from  the  actuator,  in  qui¬ 
escent  flow  around  the  hump  geometry.  This  particular  boundary  condition  was  then  used 
for  the  control  of  two  and  three-dimensional  turbulent  flow  over  the  hump.  The  actuation 
frequency  for  these  cases  were  chosen  based  on  the  physical  limitation  of  the  graphene-based 
thermoacoustic  actuator,  which  is  Helmholtz  number  ka  ~  3. 

It  was  observed  that  high-intensity  high-frequency  actuation  modifies  the  global  behav¬ 
ior  of  the  two  and  three-dimensional  flow  over  the  hump  to  achieve  drag  reduction.  In  both 
flows,  the  recirculation  zone  downstream  of  the  hump  was  elongated  in  the  streamwise  di¬ 
rection.  The  emitted  acoustic  waves  introduced  pressure  perturbations  to  the  separated  flow 
that  suppressed  the  formation  of  large-scale  shedding  vortices  by  triggering  the  formation 
of  smaller-scale  vortices,  thus  allowing  for  the  shear  layer  to  stably  convect  further  down¬ 
stream.  This  resulted  in  the  controlled  wake  having  an  elongated  recirculation  zone  with  the 
associated  low-pressure  region  shifted  downstream  leading  to  reduced  drag  on  the  hump.  We 
noted  that  this  control  mechanism  is  different  from  those  used  by  synthetic  jets  or  plasma 
actuators,  which  enhance  mixing  between  the  high  and  low-momentum  fluids  to  reduce  the 
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size  of  the  recirculation  zone  for  lower  drag.  The  present  acoustic  control  achieved  approxi¬ 
mately  10%  and  6%  of  drag  reduction  for  two  and  three-dimensional  flow,  respectively,  with 
HOdB  of  actuator  input. 

The  results  from  this  9-month  STIR  effort  appear  promising  to  aid  the  development  of 
thermoacoustic  actuator  that  can  be  locally  mounted  on  a  surface  of  arbitrary  geometry. 
The  added  attractiveness  of  this  actuator  is  the  reduced  power  input  to  deliver  the  high- 
intensity  waves  as  the  output.  The  required  power  input  is  inversely  proportional  to  the 
material  heat  capacity,  which  is  extremely  low  for  carbon-based  materials  (e.g.,  graphene  or 
carbon  nanotube  based  membranes).  To  gain  further  understanding  of  the  physics  behind 
how  locally  introduced  high-intensity  acoustic  waves  can  be  used  for  modifying  the  behavior 
of  turbulence,  additional  large-scale  simulations  and  stability  analysis  should  be  performed. 
It  would  also  be  essential  to  conduct  some  parametric  analysis  to  examine  the  influence  of 
sound  pressure  level  and  actuator  frequency  to  design  effective  and  efficient  acoustic  flow 
control  strategy.  The  results  reported  here  should  be  able  to  serve  as  a  basis  to  advance 
thermoacoustic  actuator  and  flow  control  technologies. 
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